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Abstract 

We study the Magnus effect: deflection of the trajectory of a spin- 
ning body moving in a medium. The body is rough; that is, there 
are small cavities on its surface. We concentrate on the extreme case 
of rare medium, where mutual interaction of the medium particles is 
neglected and reflections of the particles from the body's surface are 
\^J I elastic. Now the effect is due solely to multiple reflections of particles 

QQ ' in the cavities and can be studied by means of billiard theory. 

^O . We limit ourselves to the (two-dimensional) case of rough disk mov- 

z^ , ing on the plane, and to a zero-temperature medium. We calculate 

• ' the force of resistance acting on the body, as well as the moment of the 

Q ■ force, and examine how these values depend on the "kind of roughness" 

OO ! (that is, on the shape of cavities). We show that a nonzero transversal 

component of the resistance force generally appears, resulting in de- 
flection of the disk trajectory. The trajectory is determined for several 
^ , particular simple cases. 

H ' We compare our results with the previous works on Magnus effect 

in highly rarefied gases, which derive the effect from non- elastic inter- 
action of gas particles with the body. We propose another mechanism 
of creating the transversal force due to multiple collisions of particles 
with the body. 
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1 Introduction 

1.1. We are concerned with a spinning solid body moving in a homogeneous 
medium. The medium is extremely rarefied, so that the free path length of 
the medium particles is much larger than the body's size. In such a case, the 
interaction of the body with the medium can be described in terms of free 
molecular flow, where a flow of point particles falls on the body's surface; 
each particle interacts with the body but does not with other particles. There 
is no gravitation force. We suppose that the particles of the medium stay at 
rest, that is, the absolute temperature of the medium equals zero. In a frame 
of reference moving forward together with the body, we have a parallel flow 
of particles falling on the body at rest. 

Neglecting the angular momentum of particles, each particle is identified 
with a mass point that approaches the body, makes several (maybe none) 
collisions with its surface, and then goes away. All reflections are supposed 
to be perfectly elastic. 

In this paper, we restrict ourselves to the two-dimensional case. Consider 
a rough circle, that is, a set obtained from a circle by making infinitely small 
dimples on its boundary. More precisely, consider a sequence of sets Bm, 
m = 3, 4, 5, . . ., inscribed in the circle Bj.{0) of radius r centered at a point 
O. Each set B^ is invariant under the rotation by the angle 27i/m, and 
the intersection of Bm with the (27r/m)-sector AmOCm formed by two radii 
OAm and OCm is a set bounded by these radii and by a piecewise smooth 
non-self-intersecting curve contained in the triangle AmOCm and joining the 
points Am and Cm- These curves are similar for all m. A rough circle (or 
rough disk) B, is associated with such a sequence of sets B^- The curve is 
called the shape of roughness. All the values related to resistance or dynamics 
of the rough disk that are calculated below are understood as limits of the 
sequences of the corresponding values for Bm as m ^ oo. 

The roughness introduced here is uniform: it is identical at each point of 
the circle boundary. 



The rough circle is rigid and can only undergo translations or rotations. 
Denote by ip{t) the rotation angle at the time t, by uj{t) the angular velocity 
of the body, uj(t) = dip/dt, and let v(t) be the velocity of the body's center 
of mass. Let us agree to measure the rotation angle and the angular velocity 
count erclo ckwise . 
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Figure 1: A rotating rough disk in a parallel flow of particles. 



Here we consider two problems as follows: (A) determine the force of the 
medium resistance acting on the body, find the moment of this force with 
respect to the body's center of mass, and investigate their dependence on the 
shape of roughness, and (B) analyze the motion of the body in the medium, 
that is, study the functions uj{t) and v{t). Problem (A) is primary with 
respect to problem (B). In the paper, we will devote the main attention to 
problem (A), having just touched upon problem (B); we will restrict ourselves 
to deducing equations of motion and solving these equations for some simple 
particular cases. 



With each set B^ we associate a distribution of mass inside B^ such 
that the total mass M is constant and the center of mass coincides with 
the center O of the set. We also assume that the moment of inertia Im of 
Bm about O converges to a positive value J as m — *> oo. One always has 
/ < Mr'^. Denote j3 = Mr'^/I; that is, (3 is the inverse relative moment of 
inertia; we have 1 < (3 < +oo. In what follows we will pay special attention 
to the two particular cases: (a) /? = 1; that is, the mass of the rough circle 
is concentrated near its boundary, and (b) /3 = 2; the mass is distributed 
uniformly in the circle. 

The resistance force Rm{Bm,<^,uj,v) acting on the set Bm and the mo- 
ment of this force Ri^m{Bm,^p,uj,v) depend on the shape of roughness, the 
rotation angle (f, the angular velocity uj, and the velocity of translation v. 
The equations of dynamics are 

dv ^ 

M— = Rm{B.m,^,UJ,V), 

dip 

—- = UJ. 

dt 

Taking the limit m — > oo, one gets the resistance force R{B,uj,v) and the 
moment of the force Ri{B,ui,v) of the rough disk, which do not depend on 
if anymore, and the equations for the disk dynamics take the form 

M^ = R{B,uj,v), (1) 

I^ = Rj{B,co,v). (2) 



1.2. We shall see below that, generally speaking, R = R{B,u,v) is not 
collinear to v. 

If a transversal force acts on a spinning body moving in a medium, re- 
sulting in deflection of the body's trajectory, then we encounter the so called 
Magnus effect. Well-known examples of this phenomenon are dynamics of a 
golf ball and a football. The most intensively studied shapes are sphere and 
cylinder (in the second case, the cylinder is supposed to move in a direction 



orthogonal to the symmetry axis). If the direction of the transversal force 
coincides with the instantaneous velocity of the front point of the body, then 
a proper Magnus effect takes place. If these directions are opposite, then a 
reverse Magnus effect occurs. See Fig. 2 (a), (b). 
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Figure 2: (a) Magnus effect. 



(b) Reverse Magnus effect. 



There is a vast literature devoted to Magnus effect, motivated by sports 
and technology applications (see, e.g., [9],[in],Il])- Under normal conditions, 
the direct Magnus effect was found to take place. However, in certain cir- 
cumstances (for example, a well polished ball) the reverse Magnus effect can 
occur. 

Nowadays, the growing interest to aerodynamics of artificial satellites 
stimulates the study of this phenomenon in rarefied media [2j,[3j,[l]J,^12j. In 
these papers the medium is supposed to be so rarefied that a description in 
terms of free molecular flow is possible. They suggest that in this case the 
reverse Magnus effect occurs, and study this phenomenon. Moreover, in a 
recent paper the effect of Magnus force on a spinning artificial satellite is 
investigated \1]. 

In our opinion, the reverse Magnus effect in highly rarefied media is caused 
by two factors: 

(i) Non-elastic interaction of particles with the body. A part of the tan- 
gential component of the momentum of particles is transmitted to the body. 



resulting in creation of a transversal force. 

(ii) Multiple collisions of particles with the body due to the fact that the 
body's surface is not convex but contains microscopic cavities. 

In the papers [2] , |3] , [H] , |12] the impact of factor (i) is studied. Moreover, 
in these works the body is supposed to be convex and therefore factor (ii) is 
excluded. In the present paper, in contrast, we concentrate on the study of 
factor (ii). We suppose that all collisions are perfectly elastic and so exclude 
the influence of factor (i). 

We believe that the methods developed here can be extended to the three- 
dimensional case, and to the case of a finite temperature medium. Later on 
we shall see that although both reverse and direct Magnus effects are possible 
in our model, the reverse effect is the predominant one. 

Note that the limiting case of slow rotation was considered in [5], [7], [8]. 
Actually, the class of bodies studied there is rather large and includes both 
rough disks and bounded simply connected sets. It was found that the mean 
resistance force is parallel to the direction of the body's motion, and therefore 
the Magnus effect does not appear. 

In the next section, to each rough circle B we assign a measure ub charac- 
terizing the law of billiard scattering on B. The values R and Ri are defined 
to be functions of the values ub, uj, and v. In section 3 we calculate R and 
Rj for some special values of ub (and thus for some special kinds of rough 
bodies). Further, we define the set of all possible values of R, when uj and 
V are fixed and ub takes all admissible values. In other words, we answer 
the following question: what is the range of values for the resistance force 
acting on a rough unit circle? While we look through all possible shapes 
of roughness, the values of vector R cover a fixed convex two-dimensional 
set. The problem of determination of this set is formulated in terms of a 
special vector-valued Monge-Kantorovich problem and is solved numerically 
for some fixed values of the parameter A = ur/v. In section 4 we deduce 
the equations of dynamics in a convenient form and solve these equations 
for some simple particular cases. Finally, in section 5 a comparison of our 
results with the previous works on reverse Magnus effect in rarefied media is 
given. 



2 Law of billiard scattering and resistance 

2.1 Billiard scattering on a nonconvex set 

Let us first define the measure ud characterizing bilhard scattering on a 
bounded simply connected set D G M."^ with piecewise smooth boundary. 
The set d{convD) \ dD = Uj>iJj is the union of a finite or countable family 
of connected components Jj, i = 1,2,.... Each component /, is an open 
interval. (On Fig. 3, the intervals li, i = 1, 2, 3 are shown dashed.) Denote 
Jo = d{convD) ndD; in other words, Iq is the "convex part" of the boundary 
dD. Thus, dD is the disjoint union dD = Ui>o/j. 




Figure 3: Cavities on the boundary of a nonconvex set. 

Further, the set convD\D is the union of a finite or countable collection of 
its connected components. For any Jj there exists a set Qi from this collection 
such that Jj C dfli] see Fig. 3. The pair {Qi, li) will be called a cavity, and 
the interval /», the opening of the cavity. 

Denote by n^ the outer unit normal to d{convD) at the point C, € d{convD). 
On the set Jj x [— vr/2, n/2] with the coordinates (^,v^), define the measure 
fii according to the formula dfii = cos ip d^ dip, where dC, and dip stand for 
the one- dimensional Lebesgue measure. Consider the billiard in M^ \ D. Fix 
C, & li and ip G [— 7r/2, vr/2], and consider the billiard particle that starts 



to move from the point ^ with the velocity forming the angle ip with — n^. 
The particle makes one or several reflections at points of dVti \ li and then 
intersects U once again at the point ^"'" = iti.iiV)i the velocity at the mo- 
ment of intersection forming the angle ip'^ = ipl{C,,'p) with the vector n^+, 
ip^ G [— vr/2, 7r/2]. Notice that it is always vahd n^ = n^+. On Fig. 3, we 
have ip < and p^ = pf{C,, p>) > 0. 

Thus, for each i we have defined the mapping {$,,p) ^-^ i^t^vt) from a 
full measure subset of /» x [—7i/2, 7r/2] onto itself. This mapping is a bijection 
and involution, and preserves the measure /ij. In the particular case, where 
i = 0, it holds ^Q = C, and ip^ = —p. Denote by /j = | Jj| the length of /«, and 
by / = X]i>o ^i ~ \d{convD)\ the perimeter of convZ^. Introduce the notation 
D := [—7i/2, 7i/2]x [— 7r/2, it/2] and define the measures z/|^, i = 0, 1, 2, . . . 
on n as follows: i'}j{A) := j-fJ^i{{{C,,<p) '■ {p,pl{$,,'p)) G A}) for any Borel 
set A G n. In particular, the measure u^ =: uq is supported on the diagonal 
p^ = —p, and its density equals dh'Q{p,p^) = cosy? ■ 6{p + p~^). Finally, 
define ud := \ Y.i>o ^i^h- 

In a less formal way, the measure vd can be interpreted as follows. Place 
the body D in a kind of "ether" and keep it motionless. Here by ether 
we mean a homogeneous isotropic medium composed of (mutually noninter- 
acting) particles moving freely with unit velocity in all possible directions. 
When colliding with the body, the particles reflect elastically from its bound- 
ary. Thus, the particles of the ether behave like billiard particles in M^ \ D. 
Homogeneity and isotropy mean that the ether density is constant and the 
velocities of the particles are uniformly distributed in S^. 

For each particle that has reflected from D, fix the pair ((p, (/?+) consisting 
of the angle of incidence p and the angle of reflection p^ . The angle p is 
formed by the initial velocity of the particle and the vector — n^, whereas the 
angle p^ is formed by its final velocity and the vector n^+. Here ^ and ^^ are 
the points of the first and second intersection of the particle's trajectory with 
d{corwD). The distribution of the set of pairs {p,p~^) for all particles that 
have collided with D during a unit time period is described by the measure 

I'D- 

Recall that given a Borel mapping ti : X ^ Y oi two sets X C M*^^, 
Y C M'^^, with the set X being equipped with a Borel measure /x, the so- 
called push-forward measure 7r*/i on Y is defined by 7r#/i(74) := ^{tt~^{A)) 
for any Borel set A d Y . Denote by 7T^p and 7r<p+ : D -^ [— 7r/2, 7r/2] the 
projections of the square D on its horizontal and vertical sides, respectively; 



7r<^(y9, y9+) = y9, Ti^p+{ip, cp^) = (p~^. The push-forward measures 7r*z/ and vr*+z/ 
are called marginal measures for u. They are defined on [—tt/2, vr/2], and 
for any Borel set A C [— 7r/2, 7r/2] holds 

7r*iy{A) = u{A x [-7r/2, 7r/2]), 7rJ+z/(A) = z/([-7r/2, 7r/2] x A). 

Define the transformation of the square tt^ : D -^ D exchanging the coor- 
dinates if and 99+; that is, nd{^,^^) = {(p^,(p). The push-forward measure 
nfu satisfies the condition n*u{A) = ^{nci^A)) for any Borel set A C O. In 
other words, the measures irju and u are mutually symmetric with respect 
to the diagonal ip = ip~^. 

Finally, define the measure 7 on [—tt/2, 7r/2] by d'-f = cos (f dip, and denote 
by n*^^™™ the set of measures i/ on D satisfying the conditions 

7r#z/ = 7 = vrj+i/ (3) 

and 

7ifiy = u. (4) 

In other words, a generic measure from n^^^™™ is symmetric with respect to 
the diagonal (p = (p~^, and both of its marginal measures coincide with 7. One 
has z/|) G If^y™™; this can be easily deduced from the measure preserving and 
involutive properties of the mapping (^, v^) 1— > {Cty'^t)j ^^^ details see [5]. 
Hence u^ G If^^™™. The following reverse statement also holds true. 

Theorem 1. Whatever the two sets Ki C K2 C M^ such that dist(9ii'i, dK2) > 
0, the set of measures {vn : Ki G D G K2} is everywhere dense in n^^^™™ 
in the weak topology. 

The proof of this theorem can be found in [7] . 

Now consider the rough circle B generated by a sequence of sets B^- AH 
the cavities of all the sets Bm are similar, therefore ub^ does not depend 
on m and we can set by definition ub '■= ^B,n- We will see later that the 
resistance of B can be written down as a functional of i>b- The following 
theorem is obtained by a slight modification of the proof of Theorem [TJ 

Theorem 2. Whatever r > 0, the set {ub : B is a rough circle of radius r} 
is everywhere dense in If^^™™ in the weak topology. 



as columns; for instance, a vector x will be denoted by 



or (xi; X2) 



2.2 Resistance of a rough disk 

Denote v = \v\ and choose the (non-inertial) frame of reference OxiX2 in such 
a way that the direction of the axis 0x2 coincides with the direction of the 
disk motion and the origin O coincides with the disk center. In this frame of 
reference the disk stays at rest, and the flow of particles falls down on it at 
the velocity —Vq = (0; —v)"^. Here and in what follows, we represent vectors 

Xi 

X2 

Let us calculate the force R of the medium resistance and the moment of 
this force Rj with respect to O. To that end, first we consider the prelimit 
body Bm- Parameterize the opening of each cavity by the variable C, varying 
from to 1 (recall that all the cavities are identical). Denote by p the flow 
density, by (/?, the rotation angle of the cavity (that is, the external normal at 
the cavity opening equals n^ = (— sin ip; cos v^)^), and by vt. (^, ip), the final 
velocity of the particle entering the cavity at the point ^ with the velocity 
—vq- Note that At = 27i/{ujm) is the minimal time period between two 
identical positions of the rotating set B^- Then the momentum imparted to 
Bm by the particles of the flow during the time interval At equals 

2rpvAt / {-vo-v^ J ^,yD))- cos <pd^d^, (5) 

Jo J~Tv/2 ^ 

Consider the frame of reference OxiX2 having the center at the midpoint 
of the cavity opening J, the axis Oxi being parallel to /, and 0x2, codi- 
rectional with n^. That is, the frame of reference rotates jointly with the 
segment /. The change of variables from x = {xi, X2Y to x = (xi; X2Y and 
the reverse one are given by 



X = A_(^tX — r cos(7r/m) 6,^/2 and x = A^^x + r cos(7r/m) e. 



'7r/2+ajt; 



, . , COS — sin ■ , 

where Aa, = \ . , , and 



sm cos ( 



COS0 

sine/) 



Suppose now that x(t) and x(t) are the coordinates of a moving point 
in the initial and rotating frames of reference, respectively, and let v = 
{vi, V2) = dx/dt and v = [vi] V2) = dx/dt. Then 

V = A_^tV-ujA^/2-^tX and v = A^tV + ujA^/2+^tX-ujr cos^Tr/m) e^f (6) 

10 




ujr 



Figure 4: A particle falling on a cavity. 



We apply formulas ([H]) to the velocity of the particle at the two moments 
of its intersection with /. At the first moment, it holds out = ip and x = 
^e^/2+(^ + 0(1) as m — >• 00. (Here and in what follows, the estimates o(l) 
are not necessarily uniform with respect to ^ and (p.) Then the incidence 
velocity —vq takes the form 

— Vq = V {X — simp; —cosy?) + o{l) = —vg{—smx; cosx) +o(l), (7) 

where A = ur/v, 

Q = g((p, A) = V A^ — 2Asin(y9 + 1 and x = x{(p, A) 



X — simp 
arcsm — -, — r^. (8) 



q{p,X) 

As m ^ 00, the time spent by the particle in the cavity tends to zero, 
therefore the rotating frame of reference can be considered "approximately 
inertial" during that time, the velocity at the second point of intersection 
given by 

v^ =vg{- sin y; cosy f + o{l), where y = y{^,p, X) = p^{^,x{ip, X)). 

(Here ^~^{^,p), f~^{^,f) denote the mapping generated by the cavity; see 
subsection 2.1.) Applying the second formula in ([6]) and taking into account 
that X = 0(1) and ut = p + o(l), we find the velocity in the initial frame of 
reference 

^ = ^)(^'V5'-^) = vgA^{-siny; cosy)^ -vXe^ + o{l) = t;+(^, v?. A) + o(l), 

11 



where 



w^(^,V5, A) = v 



— gsm{ip + y) — X cos ip 
Q cos(v9 + y) — A sin (/? 



(9) 



Letting m — > oo in the formula ((3) for the imparted momentum and 
dividing it by At, we get the following formula for the force of resistance 
acting on the disk 



R 



Rt 
Rl 



rpv / {-Vo-v'^{^,ip,X)) COS ipd^d(p. (10) 

^0 J-tt/2 



The angular momentum transmitted to Bm by an individual particle 
equals rf^(sinx + siny) + o(l) times the mass of the particle. Summing 
the angular momenta up over all incident particles and passing to the limit 
m — >^ oo, one finds the moment of the resistance force acting on the disk, 



Ri 



r^pv 



-1 /.7r/2 



/ / vg{ip, X) {smx{ip, X) + smy{^,(p^ X)) cosipd^dif. (11) 

Jo J -it/2 



Below we will see that resistance and moment of resistance can be written 
down in the form 

(12) 



R= -rpv"^ ■ R[iyB,X], 
Ri = -r^pv^-Rj[uB,X]. 



(13) 



The dimensionless values ^[z/. A] = {Rt[i^, A]; Rl[^, X])'^ and -R/[z^, A] will be 
determined for several values of z/ and A. 

Let us consider separately the following three cases: A = l, 0<A<1, 
and A > 1, where the details of calculation are different. 



2.2.1 A 



In this case we have x = x{(p, 1) = arcsin a/(1 — smip)/2 = 7r/4— (^9/2, and so, 
the function ip ^-^ x{ip, 1) is a bijection between the intervals [— 7r/2, 7r/2] and 
[0, vr/2]. Further, one has g = g{(p,l) = a/2(1 — sinip) = 2sinx, cos(p = 
sin2x, and we obtain from (19]) 



-2 sin X cos(2x ~ y) — sin 2x 
2 sin X sin(2a; — y) — cos 2x 
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-Vo - v^ 



2v sin X 



wherefrom 

cos(2a; — y) + cos a; 

— sin(2x — y) — sin a; 

Making the change of variables {^, </?} — ^ {^i^} iii the integral in fITU]) and 
using flT^ . one gets 



R[yB, 1] = 3 



1 i-k/2 




SIV? X 



JO 



cos(2x — y) + COS a; 
— sin(2a; — y) — sin a; 



cos xd^dx. 



In this integral, y is the function of ^ and x, y = ip^{^,x). Changing 
the variables once again, {^,a;} -^ {x,y}, and taking into account that 
cos xdC,dx = di>B{x,y), we obtain 



Rli'B, 1] = 3 / / sin X 

cm 



cos(2x — y) + cosx 
- sin(2x — y) — sin x 



duB{x,y). (14) 



Here the symbol CM stays for the rectangle x G [0, vr/2], y G [— 7r/2, 7r/2]. 
The moment of the resistance force is calculated analogously, resulting in 

/•I i-n/2 

i?/[z/B, 1] = — 3 / / sin^x (sinx + siny) cosxfi^da; = 
Jo Jo 



3 // sin^ X {sinx + siny) di'B{x,y). (15) 



a 



2.2.2 0<A<1 

The second relation in ([8]) implies that for a fixed value of A, a; = x{ip,X) 
is a monotone decreasing function of ip that varies from 7r/2 to —tt/2 as (^ 
changes from — 7r/2 to 7r/2. From formula ([7]) and the first relation in (^ we 
have 



simp = X cos^ a;— sin xVl — X^ cos^ x, cos(y9 = cosx (Asinx+Vl — A^ cos^x). 



A sin X + Vl — A^ cos' x. 



1 2 ^^o2 



Denote 



( = C{x, X) = arcsin Vl — A^cos^ x; 



(16) 



13 



then one has 

cos C = -^ COS X, X + C = 71/2 — ip, C G [arccos A, 7r/2] 
and taking into account ([9]), we get 

x-y 



-vq — iF' = v{X sin a; + sin C) • 2 cos ■ 



cos (C + ^) 
in(C + ^)J 



sm 



cos if 



cosx 



A sin X + sin C = ^, 



dip d( 

dx dx 



A sin X + sin ( 
sin^ 



Using the obtained formulas, making the change of variables {^, (f} 
{^,x} in the integral ( ITO|) . and taking into account ( TT2l) . one gets 



^[^B,A] = - 



3 /'Y''/^ (Asinx + sinC)^ x-y 



OJ-tt/2 



sin^ 



cos ■ 



cos (C + ^ 
-sin(C + ^)J 



cos X d^ dx. 



Finally, the change of variables {^, x} —>■ {x, y} results in 



RWb,X] = ^ 



3 ff{\smx + sm(y x — y 



D 



sinC 



■ cos ■ 



COS (C + ^) 

-sin(C + ^). 



duB{x,y). (17) 



Recall that the symbol D denotes the square [— 7r/2, 7r/2] x [— 7r/2, 7r/2] and 

c = ax,x). 

In a similar way, from (TTT!) one gets 

3 /-y-Z^ lAsinx + sinQ^ . ^. . ,,, 

Kj = — / / (smx + smy) cos xdt,dx; 

8 Jo J-n/2 sm C 

wherefrom 

nr M 3 /"/■ (Asinx + sinC)^ 

^/i^i?,A = -- // ^— [smx + smy) duB{x,y). (18) 

8yy smC 

n 



Formulas flMj) and flTSl) are the particular cases of fITTI) and flT8|) for A = 1. 
This can be easily verified taking into account that C(x, 1) = |x|. 
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2.2.3 A > 1 



In this case, x = x{ip, A) ([8]) is not injection anymore. When ip varies from 
— 7r/2 to ipo = V^o(A) := arcsin i, the value of x monotonically decreases 
from tt/2 to Xq = Xq{X) := arccos^-, and when ip varies from ipQ to 7r/2, 
X monotonically increases from xq to it/2. Denote by </?_ := (y9_(x, A) and 
ip+ := ip^{x,\) the functions inverse to x{(p,X) on the intervals [—n/2, ipo] 
and [ipo, '/r/2], respectively. Then one has 



sm 



(p± = Acos^x ± sina;Vl — A^ cos^ 



X. 



Here and in what follows the signs "+" and "— " are related to the functions 
ip+ and (p-, respectively. The values (p+, (p-, and ( = (i^^^) (jTSl) satisfy 
the relations 



7r/2 - V5+ = X - C, vr/2 - (p_ 



x + (. 



The function ( is defined for x E [xq, tt/2] and monotonically increases from 
to 7r/2, when x changes in the interval [xq, tt/2\. 
After some algebra one gets 



cos(y9-|- sin(x =F C) 



cosx 

±^ 

dx 



cosx 



dx 



= Asinx =F sin^; 
A sin X =F sin C, 



sinC 



Q± = Asinx =F sinC; 



-vo — v^ = f (A sin X =F sin ^) ■ 2 cos 



x-y 



cos (^ T C) 
-sin(2^TC)J- 



Here the shorthand notation g± = g{ip±{x, A), A), n± = v'^iC,, 'P±{x, A), A), 
2/ = Z/(^, '/'ila;, A), A) = ¥?+(^, x) is used. ^ 

The resistance force takes the form RI^b, A] = -R_ + R+, where 



pi pTT/2 

Jo Jxo 



1 f1T/2 

{-Vo - V^) 
Jxo COSX 




COS(p± 



dip± 
dx 



cos X d^ dx 



(Asinx =F sinO^ x — y 

: cos 

sinC 2 



cos (^ T ■ 

-sin(^TC), 



cos X d^ dx. 
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Summing the integrals above and making the change of variables, one 
obtains 



RWb, A] 



cos 



x-y 



2JJ sinC 



(A sin x + 3Asinxsin (^)cosC 



cos 



sm 



x-y 
2 
x-y 



(3A sin xsin^ + sin ^)sin^ 



sm 
cos 



x-y 

2 
x-y 



duB{x,y). (19) 



Here the symbol CM stands for the rectangle [xq, 7r/2] x [— 7r/2, 7r/2]. 
The moment of the resistance force is calculated analogously. One has 
Rii^B, A] = Ri- + Ri+, where 

3 [^ r'^ cosv9± /, d^±\ , . 

Hi± = — / / g± ±— — (smx + smy) cos x at, ax = 

cosx \ ax J 



" Jo Jxo 

3 r^ r^^ {XsinxTsinCf 
8 Jo Jxo sin C 



^sin X + sin y) cos x dC, dx. 



Therefore 



„r M 3 /"V^^^ A3sin3a; + 3Asina;sin2C ^ . 

Rj[i'b,a\ = —- I I : — (sm X + sm y) cos a; (i^ ax. 

sm c, 



0^X0 



Making the change of variables, we have 

3 /■ /■ A^ sin^ X + 3A sin X sin^ ^ 
-K/[z/ij, AJ = -- // — {smx + smy) duB{x,y). (20) 

cm 



3 Magnus effect 



We are concerned here with calculation of the "normalized" resistance force 
-R[z/, A] and its moment -R/[z/, A] for some special cases of u and A, as well as 
determination of the two-dimensional set {-R[z^, A] : z/ G n^^^™™} for several 
values of A. 
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In the previous section, the formulas for i?[z/, A] and Ri[i', A] were derived, 



R[i^, A] 



RtWA] 
RlW,M 



c{x,y,X)dp(x,y) 



(21) 



n 



and 



^/h-^] = // ci{x,y,X)du{x,y), 
n 



with c and c/ given by the relations below fl22l) - fl28l) . Recall that A = 
ujr /v is the relative angular velocity and the symbol D denotes the square 
[— 7r/2, 7r/2] x [— vr/2, 7r/2]. In the formulas below, the value C = Ci^i^) is 
given by flT6|) . Xq = Xq{X) = arccos^, and x stands for the characteristic 
function. 



(a) If < A < 1 then 

3 (Asinx + sin^)^ 



c{x, y, A) 



cos 



x-y 



8 



ci{x,y,X) = 
and in particular, 

c{x,y, 1) = Ssin^x 



sinC 2 

3 (A sin X + sin (Y 



cos (C + ^ 



sinC 



sin(C + ^ 



(sin X + sin y) 



cos(2x — y) + cosx 
- sin(2x — y) — sin x 



Xx>oix,y), 



ci{x,y, 1) = -3sin^x(sina; + sin2/)xx>o(a;,l/)- 
In the limiting case A — > 0+ one has 



c{x,y,X) = -- 



9A 



sin(a; — y) 
1 + cos(x — y) 



+ 0(A), 



(22) 
(23) 

(24) 
(25) 

(26) 



ci{x,y) = -— sinx (sinx + siny) + O(A^). 
8 

(b) In the case A > 1 we have 



c(x, y, A) 



O COS t) l/\'^''^ oA- •2>\ /- 

^- < X sm X + 3A sm x sm C cos C 

2 sinC ^ s; s 



cos 



sm 



2 

x~y 
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— (3A^ sin^ a; sin C + sin^ () sin ( 



sin^ 
cos 2^ 



Xx>xoix,y), 



(27) 



cj{x,y,\) 



3 A^ sin^ X + 3A sin x sin^ ( 



4 



sinC 



sin X + sin y) Xx>xo ix,y). (28) 



3.1 Special cases of rough circles 

Consider several examples of rough circles. We present here only the final ex- 
pressions for the normalized forces and their moments calculated by formulas 
(!22|) - (l28|) : the calculation details are omitted. 



1. Circle. Consider the measure uq G n*^y™™ corresponding to the ordinary 
(without roughness) circle. The measure z/q is supported on the diago- 
nal y = — x; the density of the measure is duoi^, y) = cos a; ■ S{x + y). 
One has RtIvq,}^ = Ril^o, X] = and Rl[i^o, X] = — 1. Thus, as one 
could expect, the resistance does not depend on the angular velocity 
and is coUinear to the body's velocity. There is no Magnus effect in 
this case. 

2. Retroreflector. There exists a unique measure z/^ G H'^y™™ supported 
on the diagonal x = y; its density equals di'^{x,y) = cosx ■ 6{x — y). 
We believe that there is no rough circle generating this measure; but, 
according to Theorem[2], there exist rough circles approximating it; that 
is, there exists a sequence of rough circles -B„ such that z/^^ weakly 
converge to u^, (see also [6]). One has Rt[i^^,,\] = 37rA/8, -Rl[z/^,A] = 
—3/2 and -R/fz^*, A] = — 3A/2. Thus, the longitudinal component of the 
resistance force does not depend on the angular velocity A, while the 
transversal component and the moment of this force are proportional 
to A. 

3. Rectangular cavity. The cavity is rectangular, with the width being 
much smaller that the depth, (width) / {height of the rectangle) = e for 
e > small (see Fig. 5 (a)). We also suppose that the common length of 
the convex part of the boundary dB is very small; say, equal to e. Then 
^B = ^Tcct + o(l), where z/rect = {^0 + ^*)/'2 and o(l) denotes a measure 
that weakly converges to zero as £ ^ 0^. One can easily calculate that 
-Rrfj^rect, A] = 37rA/16, /^ifz/rect, A] = -1.25 and i?/[z/rcct, A] = -3A/4. 





Figure 5: (a) A rough disk with 
rectangular cavities. 



(b) A rough disk with 
triangular cavities. 



4. Triangular cavity. The cavity is a right isosceles triangle (see Fig. 
5(b)). The measure ub ='■ Vsi has the following support (looking like 
an inclined letter H): {a; + y = — 7r/2 : x G [— vr/2, 0]} U {|/ = a; : a; G 
[-7r/4, 7r/4]} U {a; + y = 7r/2 : a; G [0, 7r/2]}. The density of this mea- 
sure equals (iz/v( a;, y) = cosa;- f X[_z__z](x) 5(a; + y + |) + X[-^,2:](a;)- 

■b{x-y)^ X[f,f](a;)5(a: + y- f)j + \^viix\-[x{-ifi\{x)b{x ^ y ^ \)- 

- Xhf ,f ] (^) ^(^^ - y) + X[o,f ] (a;) (5(a; + y - f ) j . One has R\v^, 0+] = 

(0; -72)^ andi?/[z/v,0+] = 0; ^[z/y, 1] = (1/4 + 37r/16; 37r/16-2)^. 
The rest of the values are still unknown. 

5. Cavity realizing the product-measure. Consider the measure u^ with 
the density dv^{x.,y) = ^ cosx cosy dxdy. Evidently in this case z/g, G 
-Qsymm rjj-^g augles of incidence and of reflection are statistically inde- 
pendent; so to speak, at the moment when the particle leaves the cavity, 
it completely "forgets" its initial velocity. Here we have Rt[i'(^,X\ = 
(lOA + A^) 7r/80 for < A < 1; Riliy^, 1] = -3/4 - 7r/5 ^ -1.378; and 
RjIu^, A] = — 3A/4 for any A. The remaining values are unknown. 
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3.2 Vector-valued Monge-Kantorovich problem 

Consider the following problem: determine the set of all possible resistance 
forces that can act on a rough circle. The angular velocity of the circle is 
fixed, and the roughness varies. The force is scaled so that the resistance of 
the "ordinary circle" equals (0; —1)^. Thus, the problem can be reformulated 
in the following way: given A, find the two-dimensional set 

7^A ■.= {R[u,X] : z/GH^y'""^}. (29) 

It can be viewed as a restriction of the following more general problem: find 
the three-dimensional set 

{iR[u,X];Rj[u,X]): z/ G ff^--}. 

The latter problem is more important, but also more time-consuming, and 
is mainly postponed to the future. The only exception is the case A = 1; 
on Fig. 8 several "level sets" 7^l,c = {R[i^,l] ■ v e H^^™, i?/[z/, 1] = c} are 
depicted, suggesting an idea how the corresponding three-dimensional set 
looks like. In this case -R/[z/, 1] varies between —1.5 and 0, and the level sets 
are found for 21 values c = -1.5, -1.425, -1.35, . . . , -0.15, -0.075, 0. 

Note that the functional i?, defined on the set H'^y™™ by formula (12T|) . 
will not change if the integrand c is replaced with the symmetrized function 
c^y--(x,2/,A) = i(c(x,i/,A) + c(l/,x,A)): 

D 

Denote by 11 the set of measures i> on the square D that satisfy the only 
condition Q, that is, the set of measures with both marginals equal to 
7. For any z/ G n it holds // py"^"^ du = JJ c^y™"^ ciz^^y™"^, where z/^y°^°^ = 

i(z/ + TT*u) G n'^y'^'^. It follows that 

nx={[[c''''^"'{x,y,X)du{x,y): z/ G h}. (30) 

D 



The problem of finding TZ\ in ( 13011 is a vector-valued analogue of the 
Monge-Kantorovich problem. The difference consists in the fact that the 
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cost function, and therefore the functional, are vector valued. The set TZ\ is 
convex, since it is the image of the convex set 11 under a linear mapping. 

Note that, due to formula ([26D,c^y'°"^(a;, I/, 0+) = | (l + cos(x-y))(0; 1)^, 
therefore the problem of finding TZq+ amounts to minimizing and maximizing 

the integral I / / (1 + cos(x — y)) du(x, y) over all i/ G IT. This special Monge- 



8 

a 
Kantorovich problem was solved in [5]; the minimal and maximal values of 

the integral were found to be 0.9878... and 1.5. 

On Figures 6 and 7 we present numerical solutions of this problem for the 
values A = 0.1, 0.3 and 1, as well as the analytical solution for A = 0+. The 
case of larger A requires more involved calculation and therefore is postponed 
to the future. The method of solution is the following: for n equidistant 
vectors Cj, i = 1, . . . , n on S*^, we find the solution of the Monge-Kantorovich 
problem inf (i?[z/. A], Cj) =: Tj. Here (■ , ■) stands for the scalar product. This 
problem is reduced to the transport problem of Linear Programming and is 
solved numericallyo 

Next, the intersection of the half-planes (r, Cj) > r.i is built. It is a convex 
polygon approximating the required set 7l\, and the approximation accuracy 
increases as n increases. The value n = 100 was used in our calculations. 

On Fig. 6, the sets 7l\ are shown for A = 0"*", 0.1, 0.3, and 1. The set 7^o+ 
is the vertical segment {0} x [—1.5, —0.9878], T^o.i is the thin set with white 
interior, and T^o.s is the set with gray interior. The largest set is T^i. 

On Fig. 7, the same sets are shown in more detail. On Figures 7 (b)-(d), 
additionally, we present the regions corresponding to all possible kinds of 
roughness related to triangular cavities (and to cavities formed by combi- 
nations of different triangles), with the angles being multiples of 5°. These 
regions are colored red. For A = O"*", the corresponding region is the verti- 
cal interval {0} x [—1.42, —1] marked by a (slightly shifted) dashed line on 
Fig. 7(a). The points corresponding to special values of u are indicated by 
special signs: z/q is marked by a circle, u^, is marked by a dot, t'rect, by a 
square, z/y, by a triangle, and i/^, by a circumscribed cross. 

The part of the set TZx situated to the left of the vertical axis is related 
to resistance forces producing the (direct) Magnus effect. The part of 7l\ 
situated to the right of this axis, is related to resistance forces that cause the 



^ All the computational tests were performed on a PC Pentium IV, 2.0Ghz and 512 Mb 
RAM and using the optimization package Xpress-IVE, Version 1.19.00 with the modeler 
MOSEL. 
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reverse Magnus effect. We can see that the most part of the set (in the case 
A = 1, approximately 93.6% of the area) is situated to the right of the axis. 
This suggests that the reverse effect is more common phenomenon than the 
direct one. Actually, although Theorem 2 guarantees existence of an every- 
where dense subset of TZ\ generated by shapes of roughness, we never met a 
simple shape producing the direct Magnus effect (and thus corresponding to 
a point on the left of the vertical axis). 



■M -ffi^S^ b.5 0.75 1.0 
I . ... I .... I .... I 




Figure 6: The convex sets TZ\ with A = 0+, 0.1, 0.3 and 1 are shown. The set 
7^0+ is the vertical segment with the endpoints (0, —0.9878...) and (0, —1.5). 



22 



-1.0 -0.5 0.0 0.5 1.0 
I 



4*1 



-0.2 -0.1 0.0 0.1 0.2 

I I I I I I I I I I I I I I I I I I I I I I I 




(a) A = 



(b) A = 0.1 





(c) A = 0.3 



(d) A = 1 



Figure 7: The sets TZx with A = 0+, 0.1, 0.3, 1 are shown here separately. The 
values i?[z/. A], with z/ = i/q, t-'*, t'rect, i^v, ^® are indicated by the signs o, •, n, 
V, (g), respectively. On Fig. 7(a), the region generated by triangular cavities is 
marked by a (slightly shifted) vertical dashed line. It is the interval with the 
endpoints (0,-1) and (0,-1.42). On Fig. 7(b)-7(d), the regions generated 
by triangular cavities are colored red. 
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).0 0.25 0.5 0.75 1.0 

I I I I I I I I I I I I I I I I I I I I I I I I 




0.0 0.25 0.5 0.75 1.0 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




(a) 



(b) 



Figure 8: The "level sets" TZic = { all possible values of ^[z/, 1], with -R/[z/, 1] 
c} are shown for 21 values of c, from left to right: c = 0, —0.075, —0.15, 
—0.225, . . . , —1.425, —1.5. Figures 8(a) and 8(b) represent, respectively, 
"view from above" and "view from below" on these sets. 



4 Dynamics of a rough disk 

The motion of a spinning rough disk B is determined by the values Rt [vb , M > 
RlIi^b,^], and i?/[z/s. A]. For the sake of brevity, below we omit the fixed 
argument z/ and write R{\) instead of i?[z/. A]. Recall that the absolute value 
of the disk velocity is denoted hj v = \v\, and the angular velocity equals 
uj = Xv/r. Denote by 6 the angle the velocity makes with a fixed direction 
in an inertial frame of reference. 

Using flT^ and flT^ . one rewrites the equations of motion ([I]), ([2]) in the 
form 



dt 

de 
It 

d{Xv) 
dt 



3M 
8rpv 

3/ 



Rt{X), 
lRi{X). 
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Recall that {3 = Mr"^/! is the inverse relative moment of inertia. The 
special values that (3 can take are: /3 = 1, when the mass is concentrated on 
the disk boundary, and (3 = 2, when the mass is uniformly distributed inside 
the disk. In the intermediate case, when the mass is arbitrarily (generally 
speaking, non- uniformly) distributed inside the disk, it holds P > I. 

With the help of the change of variables 

dr = —^ dt (31 

3M ^ ' 

the equations above are transformed into the following relations: 

^ = (3Ri{\)-\Rl{\), (32) 

^ = vRl{X)^ (33) 

JO 

^ - -RAX). (34) 

Denote by s the path length of the disk; thus, ds/dt = v. From (13T|) one 
readily finds that s is proportional to r, s = |^ r. 



Below we solve the system of equations fl32|) - flMj) for the cases 1-3 
considered in subsection 3.1. Next, we determine the dynamics numerically 
for some kinds of triangular cavity. A more detailed study of dynamics in 
the general case will be addressed elsewhere. 

1. Circle. One has dX/dr = —A, dv/dr = —v and dO/dt = 0; therefore 
the circle moves straightforward. Solving these equations, one con- 
cludes that its center moves according to the equation x{t) = |^ ln(t — 
to)e + Xq, where to £ IR; e E S^ and Xq G M? are constants. Thus, hav- 
ing started the motion at some moment, the circle passes a half- line 
during an infinite time. This equation also indicates that the motion 
cannot be defined for all values of t G M. 



2. Retroreflector. Here the system (I32l)- (l34l) takes the form 

dX/dr = -3A(/3 - l)/2, dv/dr = -3t;/2, dO/dr = -SnX/S. (35) 

In the case (3 = 1, one evidently has A = const. The disk moves along 
a circumference of radius M/(7rrpA) in the direction opposite to the 
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angular velocity of rotation: if the disk rotates counterclockwise then 
its center moves clockwise along the circumference. The radius of the 
circumference is proportional to the disk mass and inversely propor- 
tional to the relative angular velocity. The path length is proportional 
to the logarithm of time, s{t) = -^ ln(t — tg)- 

In the case (3 > 1, we have s(t) = ^ ln(t — to), 9 = 9q + const • 
exp(-(/5 - 1)^ s), and X = ^{(3 - l){e - Oq). The path length once 
again depends logarithmically on the time, the relative angular velocity 
A converges to zero, and the direction 6 converges to a limiting value 
6*0; thus, the values A and 9 are exponentially decreasing functions of 
the path length and are inversely proportional to the {(3 — l)th degree 
of the time passed since a fixed moment. The trajectory of motion is a 
semibounded curve that approaches an asymptote as t -^ +oo. 

3. Rectangular cavity. Equations of motion (!32l) - (!Ml) in this case take 
the form 

dX/dr = -3X{(3 - 5/3)/4, dv/dr = -5i;/4, dO/dr = -37rA/16. 

Solving these equations, one obtains t = | ln(t — tg), v = vq e~^'^/^, A = 

^^g3r{5/3-/3)/4^ ^ = ^0 + 4(/-5/3) 6^^'^^/^-'^^^ ThuS, the path depcuds 

on t logarithmically, and the relative angular velocity and the rotation 
angle are proportional to {t — toY'^^^^ and to exp ( ^^^m '^)- 

li j3 < 5/3 then A and 9 tend to infinity, and the trajectory of the disk 
center is a converging spiral. In the case /3 > 5/3, A converges to zero, 
9 converges to a constant value, and the trajectory is a semibounded 
curve approaching an asymptote as t ^ +oo. In the case j3 = 5/3, A is 
constant, and the trajectory is a circumference of radius 2M / {rcr pX) . 

Finally, let us describe numerically some triangular cavities. It is helpful 
to consider the function g{X) = XRl{X)/Ri{X) and rewrite the equation (l32ll 
in the form 

^ = -Rj{X){g{X)-f3). (36) 

The function g{X) has been calculated for several kinds of cavity. On 
Fig. 9 (a), the two cases are illustrated where the cavity is an isosceles triangle 
with the angles (i) 30°, 120°, 30° and (ii) 60°, 60°, 60°. We see that g{X) 
monotonically increases in the case (i) and has three intervals of monotonicity 
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in the case (ii). In both cases, it is vahd Ri{\) < 0. This imphes, in the case 
(i), that the disk trajectory is a converging spiral, if /? < 1.5. li P > 1.5, 
the trajectory depends on the initial conditions and may take the form of a 
converging spiral or a curve approaching a straight line. 




a(A) 



(ii) 60", 50", 50* 




01234567 

(b) 



Figure 9: The functions (a) ^(A) = XRl{X)/Ri{X) and (b) a(A) = 4i?r(A)/A 
are shown for the triangular cavities with the angles (i) 30°, 120°, 30° and 
(ii) 60°, 60°, 60°. 



The disk behavior is more rich in the case (ii) of equilateral triangle. If 
1.38 < /3 < 1.49 then three kinds of asymptotic behavior may be realized, 
depending on the initial conditions: (I) the trajectory is a converging spiral; 
(II) the trajectory approaches a circumference; and (III) the trajectory ap- 
proaches a straight line (see Fig. 10). If 1.16 < /5 < 1.38, only two asymptotic 
behaviors of the types (I) and (II) are possible; if /3 > 1.49 then the possible 
behaviors are (I) and (III); and finally, ii (3 < 1.16, the asymptotic behavior 
is always (I). 
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/ 



Figure 10: Three kinds of asymptotic behavior of a rough disk with roughness 
formed by equilateral triangles and with 1.38 < f3 < 1.49: (I) converging 
spiral (solid line); (II) circumference (dashed line); (III) curve approaching a 
straight line (dotted line). 

In the case of triangular cavities, as our numerical evidence shows, the 
function g{X) monotonically increases for A sufficiently large and liniA^+oo 5'(A) 
+00. This implies that the trajectory is a converging spiral for appropriate 
initial conditions (namely, if the initial angular velocity is large enough). If, 
besides, j3 is large enough (that is, the mass of the disk is concentrated near 
the center), the trajectory may also be a curve approaching a straight line. 
If the function g admits intervals of monotone decrease (as for the case of 
equilateral triangle) then the trajectory may approach a circumference. The 
length of the disk path is always proportional to the logarithm of time. 
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5 Conclusions and comparison with the pre- 
vious works 

In the papers [2], [3], [11], and [12], the same problem as in the present paper 
is studied: determine the force acting on a spinning body moving through 
a rarefied gas. In [3], additionally, the moment of this force slowing down 
the body's rotation is calculated. The following shapes of the body have 
been considered: a sphere, a cylinder [3l [12] , convex bodies of revolution [3] , 
and right parallelepipeds of regular polygon section [12]. The interaction of 
the body with the gas particles is in general nonelastic: a fraction 1 — a^ of 
the incident particles is elastically refiected according to the rule "the angle 
of incidence is equal to the angle of refiection" , while the fraction a^ of the 
particles reaches thermal equilibrium with the body's surface, and is refiected 
as a Maxwellian [2] , [3] , [11] . In the paper [12], a somewhat different model of 
interaction is considered, where the reflected particles acquire a fraction a^ of 
tangential momentum of the rotating body. The transversal force results from 
the tangential friction and acts on the body in the direction corresponding to 
the reverse Magnus effect. It is remarkable that for different models and for 
different shapes of the body, the formula for the transversal force is basically 
the same. If the rotation axis is perpendicular to the direction of the body's 
motion then this force equals | arMgUov, where Mg is the mass of the gas 
displaced by the body, u is the angular velocity of the body, and v is its 
translation velocity. (Note that in [2j this formula appears in the limit of 
inflnite heat conductivity or zero gas temperature.) In [12] it was found that 
for parallelepipeds of regular n-gon section with n odd, transversal force 
depends on time, and the value of this force was determined. It is easy to 
calculate, however, that the time average of the force is equal to the same 
quantity \ arMgUJv. 

In the present paper, to the contrary, it is assumed that all collisions of 
particles with the body are absolutely elastic (that is, a,- = 0), and therefore 
there is no tangential friction. We restrict ourselves to the two-dimensional 
case and suppose that the body is a disk with small cavities on its bound- 
ary, or a rough disk. The gas has zero temperature. The Magnus effect is 
due to multiple reflections of particles in the cavities. According to (IT2|) . 
the transversal force equals \a{X) MgOJv, where A = ujr/v, Mg = vrr^p, 
a(A) = a(A, I/) = (16/37r)i?r[z/, A]/A, and z/ is the measure characterizing 
the shape of the cavities. Thus, here Mg also means the total mass of gas 
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particles displaced by the body. The function a depends on both v and A. 
In particular, a varies between —0.409 and 2 for A = 0.1, between —0.378 
and 2 for A = 0.3, and between —0.248 and 2 for A = 1. We conjecture 
that liniA^oo infjy q;(A, z/) = and liniA^oo sup^ q;(A, z/) = 2. The graphs of 
the function a (A) with v corresponding to triangular cavities with the angles 
(i) 30°, 120°, 30° and (ii) 60°, 60°, 60° are shown on Fig. 9(b). We see that 
this function significantly depends on the velocity of rotation A; in general, 
the variation of aiX) with v fixed can be more that twofold. 

In our opinion, the reverse Magnus effect appears due to the following two 
factors: (i) nonelastic interaction of particles with the body, and (ii) multiple 
reflections from cavities on the body's surface. In the papers (2].[3].[TT].[T2] 
the influence of the factor (i) has been studied, while we concentrate on 
the factor (ii). Our conclusion is that the influence of both factors is uni- 
directional, so they strengthen each other. Moreover, the formulas for the 
transversal force are similar; one should just substitute the function a(A, v^ 
for ar- We have seen that a(A, v) can be significantly greater than 1, while 
ar < 1. This can be just an artefact of two-dimensionality (in our case). 
In a forthcoming work we are planning to extend our consideration to the 
three-dimensional case and to consider media with non-zero temperature. 
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